Ideas:
| Plant | Flowers | Date | lon | lat | ele | Month | Year | julian |
|---|---|---|---|---|---|---|---|---|
| Glossoloma oblongicalyx | 4 | 2015-10-19 | -78.59093 | 0.130838 | 2270 | October | 2015 | 292 |
| Gasteranthus quitensis | 2 | 2016-10-17 | -78.59770 | 0.120070 | 1940 | October | 2016 | 291 |
| Kohleria affinis | 1 | 2016-12-13 | -78.59534 | 0.126746 | 2110 | December | 2016 | 348 |
| Columnea ciliata | 3 | 2014-02-27 | -78.59934 | 0.116682 | 1960 | February | 2014 | 58 |
| Columnea medicinalis | 1 | 2014-04-23 | -78.59372 | 0.128700 | 2130 | April | 2014 | 113 |
| Drymonia teuscheri | 3 | 2016-07-28 | -78.59245 | 0.129393 | 2200 | July | 2016 | 210 |
NEEDS TO BE FIXED.
As range
Equal probability of flowering at anytime.
## sink("model/threshold_baseline.jags")
## cat("
## model {
##
## for (x in 1:Nobs){
##
## #Observation of a flowering plant
## Y[x] ~ dbern(p[x])
## logit(p[x]) <- alpha[Plant[x]]
##
## #Residuals
## discrepancy[x] <- abs(Y[x] - p[x])
##
## #Assess Model Fit
## Ynew[x] ~ dbern(p[x])
## discrepancy.new[x]<-abs(Ynew[x] - p[x])
## }
##
##
## #Sum discrepancy
## fit<-sum(discrepancy)/Nobs
## fitnew<-sum(discrepancy.new)/Nobs
##
## #Prediction
##
## for(x in 1:Npreds){
## #predict value
##
## #Observation - probability of flowering
## prediction[x] ~ dbern(p_new[x])
## logit(p_new[x])<-alpha[PredPlant[x]]
##
## #predictive error
## pred_error[x] <- abs(Ypred[x] - prediction[x])
## }
##
##
## #Predictive Error
## fitpred<-sum(pred_error)/Npreds
##
## #Priors
##
## #Species level priors
##
## for (j in 1:Plants){
##
## #Intercept flowering probability
## alpha[j] ~ dnorm(0,0.386)
##
## }
##
## }
## ",fill=TRUE)
##
## sink()
## Compiling model graph
## Resolving undeclared variables
## Allocating nodes
## Graph information:
## Observed stochastic nodes: 2548
## Unobserved stochastic nodes: 3164
## Total graph size: 15848
##
## Initializing model
## sink("model/threshold_attraction.jags")
## cat("
## model {
##
## for (x in 1:Nobs){
##
## #Observation of a flowering plant
## Y[x] ~ dbern(p[x])
## logit(p[x]) <- alpha[Plant[x]] + e[Plant[x],Site[x]]
##
## #Residuals
## discrepancy[x] <- abs(Y[x] - p[x])
##
## #Assess Model Fit
## Ynew[x] ~ dbern(p[x])
## discrepancy.new[x]<-abs(Ynew[x] - p[x])
## }
##
## #Sum discrepancy
## fit<-sum(discrepancy)/Nobs
## fitnew<-sum(discrepancy.new)/Nobs
##
## #Prediction
##
## for(x in 1:Npreds){
## #predict value
##
## #Observation - probability of flowering
## prediction[x] ~ dbern(p_new[x])
## logit(p_new[x])<- alpha[NewPlant[x]] + e[NewPlant[x],NewSite[x]]
##
## #predictive error
## pred_error[x] <- abs(Ypred[x] - prediction[x])
## }
##
## #Predictive Error
## fitpred<-sum(pred_error)/Npreds
##
## #########################
## #autocorrelation in error
## #########################
##
## #For each of observation
## for(y in 1:Sites){
## e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
## }
##
## ##covariance among similiar species
## for(i in 1:Plants){
## for(j in 1:Plants){
## C[i,j] = exp(-lambda_cov * D[i,j])
## }
## }
##
## ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
## vCov = omega*C[,] + (1-omega) * I
## tauC=inverse(vCov*gamma)
##
## #Priors
##
## #Species level priors
##
## for (j in 1:Plants){
##
## #Intercept flowering probability
## alpha[j] ~ dnorm(0,0.386)
##
## }
## #Autocorrelation priors
## gamma = 1
##
## #Strength of covariance decay
## lambda_cov = 3
## omega = 1
## }
## ",fill=TRUE)
##
## sink()
## Compiling model graph
## Resolving undeclared variables
## Allocating nodes
## Graph information:
## Observed stochastic nodes: 2548
## Unobserved stochastic nodes: 3170
## Total graph size: 19893
##
## Initializing model
## sink("model/threshold_repulsion.jags")
## cat("
## model {
##
## for (x in 1:Nobs){
##
## #Observation of a flowering plant
## Y[x] ~ dbern(p[x])
## logit(p[x]) <- alpha[Plant[x]] + e[Plant[x],Site[x]]
##
## #Residuals
## discrepancy[x] <- abs(Y[x] - p[x])
##
## #Assess Model Fit
## Ynew[x] ~ dbern(p[x])
## discrepancy.new[x]<-abs(Ynew[x] - p[x])
## }
##
##
## #Sum discrepancy
## fit<-sum(discrepancy)/Nobs
## fitnew<-sum(discrepancy.new)/Nobs
##
## #Prediction
##
## for(x in 1:Npreds){
## #predict value
##
## #Observation - probability of flowering
## prediction[x] ~ dbern(p_new[x])
## logit(p_new[x])<- alpha[NewPlant[x]] + e[NewPlant[x],NewSite[x]]
##
## #predictive error
## pred_error[x] <- abs(Ypred[x] - prediction[x])
## }
##
## #Predictive Error
## fitpred<-sum(pred_error)/Npreds
##
## #########################
## #autocorrelation in error
## #########################
##
## #For each of observation
## for(y in 1:Sites){
## e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
## }
##
## ##covariance among similiar species
## for(i in 1:Plants){
## for(j in 1:Plants){
## C[i,j] = exp(-lambda_cov * D[i,j])
## }
## }
##
## ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
## vCov = omega*C[,] + (1-omega) * I
## tauC=vCov*gamma
##
## #Priors
##
## #Species level priors
##
## for (j in 1:Plants){
##
## #Intercept flowering probability
## alpha[j] ~ dnorm(0,0.386)
##
## }
## #Autocorrelation priors
## gamma = 5
##
## #Strength of covariance decay
## lambda_cov = 2
## omega = 1
## }
## ",fill=TRUE)
##
## sink()
## Compiling model graph
## Resolving undeclared variables
## Allocating nodes
## Graph information:
## Observed stochastic nodes: 2548
## Unobserved stochastic nodes: 3170
## Total graph size: 19893
##
## Initializing model
## sink("model/threshold_attraction.jags")
## cat("
## model {
##
## for (x in 1:Nobs){
##
## #Observation of a flowering plant
## Y[x] ~ dbern(p[x])
## logit(p[x]) <- alpha[Plant[x]] + e[Plant[x],Site[x]]
##
## #Residuals
## discrepancy[x] <- abs(Y[x] - p[x])
##
## #Assess Model Fit
## Ynew[x] ~ dbern(p[x])
## discrepancy.new[x]<-abs(Ynew[x] - p[x])
## }
##
## #Sum discrepancy
## fit<-sum(discrepancy)/Nobs
## fitnew<-sum(discrepancy.new)/Nobs
##
## #Prediction
##
## for(x in 1:Npreds){
## #predict value
##
## #Observation - probability of flowering
## prediction[x] ~ dbern(p_new[x])
## logit(p_new[x])<- alpha[NewPlant[x]] + e[NewPlant[x],NewSite[x]]
##
## #predictive error
## pred_error[x] <- abs(Ypred[x] - prediction[x])
## }
##
## #Predictive Error
## fitpred<-sum(pred_error)/Npreds
##
## #########################
## #autocorrelation in error
## #########################
##
## #For each of observation
## for(y in 1:Sites){
## e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
## }
##
## ##covariance among similiar species
## for(i in 1:Plants){
## for(j in 1:Plants){
## C[i,j] = exp(-lambda_cov * D[i,j])
## }
## }
##
## ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
## vCov = omega*C[,] + (1-omega) * I
## tauC=inverse(vCov*gamma)
##
## #Priors
##
## #Species level priors
##
## for (j in 1:Plants){
##
## #Intercept flowering probability
## alpha[j] ~ dnorm(0,0.386)
##
## }
## #Autocorrelation priors
## gamma = 1
##
## #Strength of covariance decay
## lambda_cov = 3
## omega = 1
## }
## ",fill=TRUE)
##
## sink()
## Compiling model graph
## Resolving undeclared variables
## Allocating nodes
## Graph information:
## Observed stochastic nodes: 2548
## Unobserved stochastic nodes: 3170
## Total graph size: 20049
##
## Initializing model
## sink("model/threshold_repulsion.jags")
## cat("
## model {
##
## for (x in 1:Nobs){
##
## #Observation of a flowering plant
## Y[x] ~ dbern(p[x])
## logit(p[x]) <- alpha[Plant[x]] + e[Plant[x],Site[x]]
##
## #Residuals
## discrepancy[x] <- abs(Y[x] - p[x])
##
## #Assess Model Fit
## Ynew[x] ~ dbern(p[x])
## discrepancy.new[x]<-abs(Ynew[x] - p[x])
## }
##
##
## #Sum discrepancy
## fit<-sum(discrepancy)/Nobs
## fitnew<-sum(discrepancy.new)/Nobs
##
## #Prediction
##
## for(x in 1:Npreds){
## #predict value
##
## #Observation - probability of flowering
## prediction[x] ~ dbern(p_new[x])
## logit(p_new[x])<- alpha[NewPlant[x]] + e[NewPlant[x],NewSite[x]]
##
## #predictive error
## pred_error[x] <- abs(Ypred[x] - prediction[x])
## }
##
## #Predictive Error
## fitpred<-sum(pred_error)/Npreds
##
## #########################
## #autocorrelation in error
## #########################
##
## #For each of observation
## for(y in 1:Sites){
## e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
## }
##
## ##covariance among similiar species
## for(i in 1:Plants){
## for(j in 1:Plants){
## C[i,j] = exp(-lambda_cov * D[i,j])
## }
## }
##
## ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
## vCov = omega*C[,] + (1-omega) * I
## tauC=vCov*gamma
##
## #Priors
##
## #Species level priors
##
## for (j in 1:Plants){
##
## #Intercept flowering probability
## alpha[j] ~ dnorm(0,0.386)
##
## }
## #Autocorrelation priors
## gamma = 5
##
## #Strength of covariance decay
## lambda_cov = 2
## omega = 1
## }
## ",fill=TRUE)
##
## sink()
## Compiling model graph
## Resolving undeclared variables
## Allocating nodes
## Graph information:
## Observed stochastic nodes: 2548
## Unobserved stochastic nodes: 3170
## Total graph size: 20049
##
## Initializing model
## sink("model/threshold_attraction.jags")
## cat("
## model {
##
## for (x in 1:Nobs){
##
## #Observation of a flowering plant
## Y[x] ~ dbern(p[x])
## logit(p[x]) <- alpha[Plant[x]] + e[Plant[x],Site[x]]
##
## #Residuals
## discrepancy[x] <- abs(Y[x] - p[x])
##
## #Assess Model Fit
## Ynew[x] ~ dbern(p[x])
## discrepancy.new[x]<-abs(Ynew[x] - p[x])
## }
##
## #Sum discrepancy
## fit<-sum(discrepancy)/Nobs
## fitnew<-sum(discrepancy.new)/Nobs
##
## #Prediction
##
## for(x in 1:Npreds){
## #predict value
##
## #Observation - probability of flowering
## prediction[x] ~ dbern(p_new[x])
## logit(p_new[x])<- alpha[NewPlant[x]] + e[NewPlant[x],NewSite[x]]
##
## #predictive error
## pred_error[x] <- abs(Ypred[x] - prediction[x])
## }
##
## #Predictive Error
## fitpred<-sum(pred_error)/Npreds
##
## #########################
## #autocorrelation in error
## #########################
##
## #For each of observation
## for(y in 1:Sites){
## e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
## }
##
## ##covariance among similiar species
## for(i in 1:Plants){
## for(j in 1:Plants){
## C[i,j] = exp(-lambda_cov * D[i,j])
## }
## }
##
## ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
## vCov = omega*C[,] + (1-omega) * I
## tauC=inverse(vCov*gamma)
##
## #Priors
##
## #Species level priors
##
## for (j in 1:Plants){
##
## #Intercept flowering probability
## alpha[j] ~ dnorm(0,0.386)
##
## }
## #Autocorrelation priors
## gamma = 1
##
## #Strength of covariance decay
## lambda_cov = 3
## omega = 1
## }
## ",fill=TRUE)
##
## sink()
## Compiling model graph
## Resolving undeclared variables
## Allocating nodes
## Graph information:
## Observed stochastic nodes: 2548
## Unobserved stochastic nodes: 3170
## Total graph size: 20049
##
## Initializing model
## sink("model/threshold_repulsion.jags")
## cat("
## model {
##
## for (x in 1:Nobs){
##
## #Observation of a flowering plant
## Y[x] ~ dbern(p[x])
## logit(p[x]) <- alpha[Plant[x]] + e[Plant[x],Site[x]]
##
## #Residuals
## discrepancy[x] <- abs(Y[x] - p[x])
##
## #Assess Model Fit
## Ynew[x] ~ dbern(p[x])
## discrepancy.new[x]<-abs(Ynew[x] - p[x])
## }
##
##
## #Sum discrepancy
## fit<-sum(discrepancy)/Nobs
## fitnew<-sum(discrepancy.new)/Nobs
##
## #Prediction
##
## for(x in 1:Npreds){
## #predict value
##
## #Observation - probability of flowering
## prediction[x] ~ dbern(p_new[x])
## logit(p_new[x])<- alpha[NewPlant[x]] + e[NewPlant[x],NewSite[x]]
##
## #predictive error
## pred_error[x] <- abs(Ypred[x] - prediction[x])
## }
##
## #Predictive Error
## fitpred<-sum(pred_error)/Npreds
##
## #########################
## #autocorrelation in error
## #########################
##
## #For each of observation
## for(y in 1:Sites){
## e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
## }
##
## ##covariance among similiar species
## for(i in 1:Plants){
## for(j in 1:Plants){
## C[i,j] = exp(-lambda_cov * D[i,j])
## }
## }
##
## ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
## vCov = omega*C[,] + (1-omega) * I
## tauC=vCov*gamma
##
## #Priors
##
## #Species level priors
##
## for (j in 1:Plants){
##
## #Intercept flowering probability
## alpha[j] ~ dnorm(0,0.386)
##
## }
## #Autocorrelation priors
## gamma = 5
##
## #Strength of covariance decay
## lambda_cov = 2
## omega = 1
## }
## ",fill=TRUE)
##
## sink()
## Compiling model graph
## Resolving undeclared variables
## Allocating nodes
## Graph information:
## Observed stochastic nodes: 2548
## Unobserved stochastic nodes: 3170
## Total graph size: 20049
##
## Initializing model
# One example
Glossoloma purpureum
## # A tibble: 7 x 2
## Model p
## <chr> <dbl>
## 1 baseline 0.528
## 2 interaction_attraction 0.974
## 3 interaction_repulsion 0.978
## 4 phylogenetic_attraction 0.994
## 5 phylogenetic_repulsion 0.98
## 6 trait_attraction 0.978
## 7 trait_repulsion 0.958
Without baseline
## # A tibble: 6 x 2
## Model p
## <chr> <dbl>
## 1 interaction_attraction 0.974
## 2 interaction_repulsion 0.978
## 3 phylogenetic_attraction 0.994
## 4 phylogenetic_repulsion 0.98
## 5 trait_attraction 0.978
## 6 trait_repulsion 0.958
| Model | mean | lower | upper |
|---|---|---|---|
| trait_attraction | 0.2059569 | 0.1995477 | 0.2119726 |
| phylogenetic_attraction | 0.2046194 | 0.1980780 | 0.2109580 |
| interaction_attraction | 0.2043841 | 0.1977721 | 0.2109840 |
| phylogenetic_repulsion | 0.1973214 | 0.1911837 | 0.2035766 |
| interaction_repulsion | 0.1971715 | 0.1904616 | 0.2031607 |
| trait_repulsion | 0.1964914 | 0.1900542 | 0.2028746 |
| Model | mean | lower | upper |
|---|---|---|---|
| baseline | 0.2434784 | 0.2192691 | 0.2674419 |
| phylogenetic_attraction | 0.2126379 | 0.1910299 | 0.2358804 |
| trait_attraction | 0.2124917 | 0.1926080 | 0.2358804 |
| interaction_attraction | 0.2110764 | 0.1877076 | 0.2342193 |
| interaction_repulsion | 0.2098140 | 0.1876246 | 0.2358804 |
| trait_repulsion | 0.2097043 | 0.1877076 | 0.2342193 |
| phylogenetic_repulsion | 0.2084718 | 0.1843854 | 0.2342193 |
| Plant | Month | Observed | baseline | interaction_attraction | interaction_repulsion | phylogenetic_attraction | phylogenetic_repulsion | trait_attraction | trait_repulsion |
|---|---|---|---|---|---|---|---|---|---|
| Besleria solanoides | Apr | 13.3 | 17.0 | 15.7 | 16.9 | 16.7 | 16.8 | 16.1 | 16.3 |
| Besleria solanoides | Aug | 5.9 | 16.6 | 16.5 | 16.5 | 16.4 | 16.5 | 15.6 | 17.7 |
| Besleria solanoides | Dec | 12.5 | 17.8 | 16.6 | 17.6 | 15.9 | 18.3 | 16.1 | 17.9 |
| Besleria solanoides | Feb | 0.0 | 16.6 | 16.0 | 16.6 | 15.5 | 16.0 | 15.8 | 16.0 |
| Besleria solanoides | Jan | 0.0 | 17.3 | 16.8 | 18.1 | 17.6 | 17.1 | 17.9 | 16.3 |
| Besleria solanoides | Jul | 5.6 | 16.8 | 14.9 | 15.3 | 15.0 | 15.8 | 14.6 | 14.8 |
| Besleria solanoides | Jun | 5.0 | 16.1 | 15.2 | 15.8 | 15.5 | 15.4 | 15.5 | 16.6 |
| Besleria solanoides | Mar | 0.0 | 17.4 | 16.8 | 17.4 | 17.5 | 17.1 | 17.4 | 17.3 |
| Besleria solanoides | May | 0.0 | 16.1 | 16.5 | 17.4 | 17.2 | 16.5 | 16.4 | 17.1 |
| Besleria solanoides | Nov | 33.3 | 16.4 | 16.1 | 17.5 | 17.3 | 15.6 | 15.9 | 17.6 |
| Besleria solanoides | Oct | 55.6 | 16.7 | 16.7 | 17.8 | 17.9 | 16.9 | 17.4 | 16.9 |
| Besleria solanoides | Sep | 50.0 | 15.5 | 16.3 | 15.8 | 15.5 | 16.4 | 15.5 | 15.7 |
| Columnea ciliata | Apr | 46.7 | 15.8 | 17.8 | 18.8 | 17.8 | 19.3 | 17.4 | 18.9 |
| Columnea ciliata | Aug | 0.0 | 16.6 | 13.5 | 12.7 | 14.2 | 13.0 | 13.9 | 12.4 |
| Columnea ciliata | Dec | 0.0 | 15.4 | 16.2 | 15.3 | 16.7 | 15.3 | 14.6 | 15.2 |
| Columnea ciliata | Feb | 43.8 | 16.2 | 15.5 | 16.5 | 16.2 | 16.2 | 15.8 | 15.6 |
| Columnea ciliata | Jan | 25.0 | 15.9 | 19.0 | 19.2 | 18.0 | 19.6 | 18.3 | 19.3 |
| Columnea ciliata | Jul | 5.6 | 16.1 | 16.4 | 17.4 | 16.8 | 17.1 | 16.7 | 16.4 |
| Columnea ciliata | Jun | 10.0 | 15.7 | 15.8 | 17.0 | 15.9 | 15.7 | 16.2 | 15.5 |
| Columnea ciliata | Mar | 66.7 | 16.1 | 17.9 | 19.4 | 17.8 | 19.9 | 18.3 | 20.2 |
| Columnea ciliata | May | 5.9 | 16.6 | 15.9 | 17.0 | 15.5 | 16.9 | 16.8 | 16.6 |
| Columnea ciliata | Nov | 0.0 | 17.1 | 16.6 | 16.4 | 15.8 | 16.6 | 16.0 | 15.5 |
| Columnea ciliata | Oct | 0.0 | 16.4 | 15.7 | 14.8 | 15.1 | 15.3 | 15.2 | 14.5 |
| Columnea ciliata | Sep | 0.0 | 15.7 | 15.6 | 16.3 | 16.2 | 16.6 | 16.1 | 15.8 |
| Columnea kucyniakii | Apr | 6.7 | 8.8 | 8.6 | 8.9 | 7.9 | 8.9 | 8.4 | 8.2 |
| Columnea kucyniakii | Aug | 0.0 | 8.8 | 8.3 | 9.3 | 8.4 | 9.4 | 8.8 | 9.2 |
| Columnea kucyniakii | Dec | 12.5 | 8.5 | 9.0 | 9.5 | 9.8 | 8.9 | 9.3 | 9.5 |
| Columnea kucyniakii | Feb | 43.8 | 8.1 | 8.2 | 9.6 | 8.4 | 8.5 | 8.4 | 9.1 |
| Columnea kucyniakii | Jan | 25.0 | 8.9 | 9.5 | 9.3 | 8.4 | 10.2 | 9.2 | 9.8 |
| Columnea kucyniakii | Jul | 0.0 | 8.7 | 7.9 | 8.8 | 7.6 | 8.3 | 8.2 | 8.4 |
| Columnea kucyniakii | Jun | 0.0 | 9.2 | 8.1 | 8.7 | 8.9 | 8.5 | 8.1 | 8.8 |
| Columnea kucyniakii | Mar | 16.7 | 8.8 | 8.4 | 10.2 | 9.0 | 9.8 | 9.0 | 8.9 |
| Columnea kucyniakii | May | 0.0 | 8.4 | 8.7 | 9.6 | 9.3 | 9.6 | 10.1 | 8.8 |
| Columnea kucyniakii | Nov | 11.1 | 8.4 | 8.5 | 8.1 | 7.8 | 8.5 | 7.6 | 9.6 |
| Columnea kucyniakii | Oct | 0.0 | 8.4 | 8.6 | 8.6 | 8.3 | 9.3 | 7.9 | 9.1 |
| Columnea kucyniakii | Sep | 0.0 | 8.8 | 8.2 | 9.6 | 8.2 | 8.9 | 9.1 | 8.9 |
| Columnea medicinalis | Apr | 0.0 | 15.7 | 14.0 | 14.8 | 14.1 | 14.3 | 14.0 | 14.6 |
| Columnea medicinalis | Aug | 23.5 | 14.1 | 16.8 | 16.0 | 15.6 | 16.6 | 16.2 | 15.6 |
| Columnea medicinalis | Dec | 0.0 | 14.6 | 16.4 | 16.2 | 15.8 | 16.1 | 16.6 | 15.9 |
| Columnea medicinalis | Feb | 12.5 | 14.0 | 14.2 | 14.5 | 14.1 | 14.1 | 14.1 | 14.2 |
| Columnea medicinalis | Jan | 16.7 | 13.8 | 13.9 | 15.0 | 15.3 | 14.7 | 14.0 | 15.5 |
| Columnea medicinalis | Jul | 27.8 | 14.3 | 13.2 | 13.0 | 13.2 | 13.0 | 12.4 | 13.2 |
| Columnea medicinalis | Jun | 20.0 | 14.6 | 14.5 | 14.3 | 14.7 | 14.9 | 14.9 | 14.6 |
| Columnea medicinalis | Mar | 8.3 | 14.3 | 14.0 | 15.2 | 14.8 | 15.4 | 13.4 | 15.0 |
| Columnea medicinalis | May | 11.8 | 14.4 | 16.2 | 16.8 | 16.3 | 15.5 | 15.2 | 16.6 |
| Columnea medicinalis | Nov | 11.1 | 14.1 | 15.5 | 14.6 | 14.6 | 15.5 | 15.1 | 14.7 |
| Columnea medicinalis | Oct | 22.2 | 13.9 | 15.8 | 15.8 | 14.5 | 16.1 | 16.1 | 15.1 |
| Columnea medicinalis | Sep | 5.0 | 13.8 | 14.3 | 14.6 | 14.5 | 14.1 | 14.0 | 14.5 |
| Columnea picta | Apr | 13.3 | 14.0 | 14.4 | 16.0 | 14.8 | 15.8 | 14.4 | 15.6 |
| Columnea picta | Aug | 5.9 | 14.3 | 12.2 | 11.3 | 12.1 | 11.5 | 11.8 | 11.4 |
| Columnea picta | Dec | 12.5 | 14.3 | 14.7 | 14.4 | 15.0 | 13.6 | 13.9 | 13.9 |
| Columnea picta | Feb | 6.2 | 14.2 | 14.0 | 15.6 | 14.5 | 14.8 | 14.6 | 13.9 |
| Columnea picta | Jan | 16.7 | 14.5 | 14.6 | 14.8 | 15.1 | 14.7 | 14.7 | 15.4 |
| Columnea picta | Jul | 16.7 | 14.0 | 14.0 | 16.1 | 15.2 | 16.0 | 14.9 | 15.2 |
| Columnea picta | Jun | 25.0 | 14.5 | 14.0 | 14.8 | 13.8 | 14.4 | 13.9 | 14.3 |
| Columnea picta | Mar | 0.0 | 14.8 | 14.5 | 14.9 | 14.2 | 15.4 | 14.0 | 14.6 |
| Columnea picta | May | 35.3 | 14.1 | 17.1 | 17.9 | 18.0 | 17.6 | 17.4 | 17.4 |
| Columnea picta | Nov | 0.0 | 13.9 | 10.6 | 10.6 | 10.5 | 10.3 | 9.1 | 10.4 |
| Columnea picta | Oct | 0.0 | 13.7 | 10.4 | 9.1 | 9.5 | 9.5 | 9.6 | 9.6 |
| Columnea picta | Sep | 20.0 | 14.0 | 15.5 | 16.2 | 15.9 | 16.8 | 15.7 | 16.5 |
| Columnea strigosa | Apr | 0.0 | 12.2 | 11.5 | 11.5 | 11.1 | 10.8 | 11.6 | 11.4 |
| Columnea strigosa | Aug | 11.8 | 12.3 | 14.6 | 15.0 | 14.2 | 15.6 | 15.3 | 15.4 |
| Columnea strigosa | Dec | 0.0 | 11.5 | 12.2 | 11.2 | 11.6 | 11.3 | 12.7 | 11.3 |
| Columnea strigosa | Feb | 12.5 | 12.6 | 12.4 | 13.9 | 12.7 | 13.2 | 13.1 | 13.7 |
| Columnea strigosa | Jan | 8.3 | 11.6 | 9.2 | 10.1 | 9.5 | 9.5 | 9.8 | 9.5 |
| Columnea strigosa | Jul | 27.8 | 12.7 | 12.8 | 14.3 | 13.9 | 13.8 | 13.3 | 14.7 |
| Columnea strigosa | Jun | 40.0 | 12.0 | 14.0 | 14.0 | 12.7 | 13.1 | 13.9 | 13.8 |
| Columnea strigosa | Mar | 0.0 | 12.0 | 8.9 | 10.1 | 9.8 | 9.3 | 10.0 | 9.2 |
| Columnea strigosa | May | 23.5 | 12.1 | 10.4 | 9.9 | 10.1 | 9.6 | 10.2 | 9.7 |
| Columnea strigosa | Nov | 0.0 | 12.3 | 13.9 | 14.9 | 13.4 | 15.3 | 14.2 | 14.3 |
| Columnea strigosa | Oct | 0.0 | 12.6 | 14.5 | 14.2 | 13.2 | 14.6 | 14.3 | 15.0 |
| Columnea strigosa | Sep | 0.0 | 11.7 | 12.4 | 12.0 | 11.9 | 11.9 | 12.4 | 12.7 |
| Drymonia collegarum | Apr | 13.3 | 12.8 | 13.2 | 12.9 | 13.7 | 13.1 | 12.9 | 13.5 |
| Drymonia collegarum | Aug | 17.6 | 13.2 | 14.9 | 14.9 | 16.1 | 14.2 | 15.6 | 13.8 |
| Drymonia collegarum | Dec | 25.0 | 14.2 | 15.1 | 15.0 | 16.0 | 16.2 | 15.6 | 15.0 |
| Drymonia collegarum | Feb | 6.2 | 13.3 | 13.1 | 12.8 | 12.5 | 12.4 | 13.2 | 12.4 |
| Drymonia collegarum | Jan | 16.7 | 13.1 | 15.2 | 15.2 | 13.9 | 15.0 | 13.4 | 15.0 |
| Drymonia collegarum | Jul | 11.1 | 13.0 | 12.0 | 11.5 | 12.5 | 10.9 | 12.0 | 11.5 |
| Drymonia collegarum | Jun | 15.0 | 13.6 | 13.2 | 13.9 | 13.3 | 13.3 | 13.2 | 12.8 |
| Drymonia collegarum | Mar | 8.3 | 13.2 | 15.8 | 15.2 | 14.1 | 15.7 | 13.2 | 15.6 |
| Drymonia collegarum | May | 23.5 | 13.5 | 15.5 | 14.7 | 13.8 | 15.2 | 13.9 | 14.8 |
| Drymonia collegarum | Nov | 11.1 | 14.3 | 13.7 | 14.4 | 14.8 | 13.4 | 14.3 | 13.5 |
| Drymonia collegarum | Oct | 5.6 | 13.5 | 13.2 | 14.7 | 14.1 | 13.9 | 14.4 | 14.0 |
| Drymonia collegarum | Sep | 10.0 | 13.9 | 13.3 | 12.9 | 13.0 | 13.1 | 12.4 | 13.0 |
| Drymonia tenuis | Apr | 13.3 | 13.9 | 14.7 | 15.1 | 14.6 | 15.5 | 14.3 | 15.4 |
| Drymonia tenuis | Aug | 11.8 | 14.1 | 11.0 | 13.7 | 13.3 | 11.8 | 13.6 | 11.6 |
| Drymonia tenuis | Dec | 0.0 | 14.1 | 13.0 | 13.3 | 13.1 | 13.1 | 13.9 | 12.8 |
| Drymonia tenuis | Feb | 18.8 | 14.1 | 14.8 | 14.0 | 14.6 | 14.9 | 13.2 | 14.5 |
| Drymonia tenuis | Jan | 16.7 | 14.8 | 15.8 | 17.2 | 17.0 | 17.0 | 15.8 | 17.2 |
| Drymonia tenuis | Jul | 16.7 | 13.9 | 13.8 | 13.8 | 13.6 | 13.9 | 12.7 | 14.7 |
| Drymonia tenuis | Jun | 20.0 | 14.6 | 13.8 | 13.0 | 13.8 | 13.7 | 13.2 | 13.8 |
| Drymonia tenuis | Mar | 25.0 | 14.2 | 17.2 | 17.0 | 17.1 | 17.0 | 15.2 | 17.5 |
| Drymonia tenuis | May | 23.5 | 13.4 | 14.6 | 13.1 | 13.9 | 12.5 | 13.5 | 13.1 |
| Drymonia tenuis | Nov | 0.0 | 14.8 | 13.6 | 16.1 | 15.5 | 14.8 | 15.6 | 16.0 |
| Drymonia tenuis | Oct | 5.6 | 14.3 | 13.9 | 16.5 | 15.3 | 15.2 | 16.0 | 14.5 |
| Drymonia tenuis | Sep | 5.0 | 13.5 | 13.3 | 13.3 | 13.4 | 13.7 | 13.3 | 13.5 |
| Drymonia teuscheri | Apr | 26.7 | 16.4 | 15.9 | 15.4 | 16.3 | 16.1 | 16.4 | 16.6 |
| Drymonia teuscheri | Aug | 23.5 | 16.5 | 19.3 | 19.7 | 19.4 | 20.1 | 18.8 | 20.7 |
| Drymonia teuscheri | Dec | 0.0 | 16.4 | 17.3 | 18.1 | 17.3 | 17.4 | 18.2 | 18.0 |
| Drymonia teuscheri | Feb | 6.2 | 16.7 | 16.5 | 16.8 | 16.8 | 16.9 | 17.2 | 17.9 |
| Drymonia teuscheri | Jan | 16.7 | 16.7 | 12.7 | 14.4 | 14.1 | 13.2 | 15.7 | 13.3 |
| Drymonia teuscheri | Jul | 22.2 | 17.5 | 17.3 | 16.3 | 15.9 | 17.8 | 16.7 | 17.6 |
| Drymonia teuscheri | Jun | 15.0 | 17.0 | 17.4 | 17.4 | 17.5 | 18.2 | 17.1 | 17.5 |
| Drymonia teuscheri | Mar | 16.7 | 16.7 | 13.2 | 14.4 | 15.1 | 13.2 | 14.1 | 13.2 |
| Drymonia teuscheri | May | 29.4 | 17.0 | 14.7 | 16.1 | 16.2 | 15.1 | 16.0 | 14.9 |
| Drymonia teuscheri | Nov | 11.1 | 16.3 | 18.8 | 18.8 | 18.1 | 18.9 | 18.1 | 20.0 |
| Drymonia teuscheri | Oct | 11.1 | 16.6 | 18.8 | 19.1 | 18.8 | 19.8 | 18.6 | 19.2 |
| Drymonia teuscheri | Sep | 10.0 | 16.6 | 16.8 | 16.4 | 16.6 | 17.1 | 16.2 | 16.8 |
| Gasteranthus lateralis | Apr | 13.3 | 14.0 | 13.3 | 14.4 | 13.5 | 13.7 | 14.8 | 13.8 |
| Gasteranthus lateralis | Aug | 17.6 | 14.0 | 16.0 | 16.0 | 15.6 | 16.4 | 15.1 | 16.8 |
| Gasteranthus lateralis | Dec | 25.0 | 14.6 | 15.2 | 15.0 | 14.5 | 14.1 | 15.3 | 15.8 |
| Gasteranthus lateralis | Feb | 0.0 | 14.5 | 14.2 | 14.2 | 14.1 | 14.3 | 13.9 | 14.1 |
| Gasteranthus lateralis | Jan | 33.3 | 14.8 | 13.0 | 12.8 | 12.9 | 12.9 | 14.1 | 12.2 |
| Gasteranthus lateralis | Jul | 0.0 | 13.6 | 13.0 | 14.6 | 13.5 | 13.9 | 13.1 | 13.8 |
| Gasteranthus lateralis | Jun | 0.0 | 14.3 | 14.2 | 15.2 | 13.9 | 14.8 | 13.9 | 14.6 |
| Gasteranthus lateralis | Mar | 8.3 | 13.7 | 13.8 | 13.3 | 12.8 | 12.7 | 13.7 | 11.8 |
| Gasteranthus lateralis | May | 5.9 | 13.3 | 13.6 | 14.1 | 14.4 | 14.0 | 13.7 | 13.0 |
| Gasteranthus lateralis | Nov | 44.4 | 14.5 | 14.8 | 14.8 | 15.0 | 16.1 | 14.6 | 15.7 |
| Gasteranthus lateralis | Oct | 33.3 | 14.0 | 14.6 | 15.8 | 14.8 | 15.4 | 15.7 | 15.4 |
| Gasteranthus lateralis | Sep | 10.0 | 14.4 | 13.4 | 14.4 | 14.3 | 14.0 | 14.0 | 14.3 |
| Gasteranthus quitensis | Apr | 0.0 | 11.7 | 12.1 | 13.0 | 12.3 | 12.4 | 11.6 | 12.2 |
| Gasteranthus quitensis | Aug | 17.6 | 12.3 | 11.2 | 10.0 | 10.3 | 10.0 | 10.4 | 11.0 |
| Gasteranthus quitensis | Dec | 12.5 | 11.9 | 11.4 | 10.9 | 10.2 | 11.0 | 11.5 | 11.5 |
| Gasteranthus quitensis | Feb | 18.8 | 12.0 | 11.4 | 13.8 | 12.3 | 12.4 | 12.7 | 12.5 |
| Gasteranthus quitensis | Jan | 33.3 | 12.3 | 12.8 | 14.4 | 14.2 | 13.1 | 14.1 | 14.1 |
| Gasteranthus quitensis | Jul | 5.6 | 11.7 | 11.8 | 12.7 | 11.8 | 12.0 | 11.7 | 12.9 |
| Gasteranthus quitensis | Jun | 0.0 | 11.4 | 11.1 | 11.8 | 11.2 | 11.3 | 10.8 | 11.5 |
| Gasteranthus quitensis | Mar | 0.0 | 11.5 | 13.0 | 14.8 | 14.0 | 13.3 | 13.8 | 14.6 |
| Gasteranthus quitensis | May | 0.0 | 12.2 | 11.6 | 11.5 | 11.6 | 11.0 | 11.3 | 12.4 |
| Gasteranthus quitensis | Nov | 11.1 | 12.5 | 12.2 | 11.5 | 12.5 | 10.8 | 11.6 | 11.1 |
| Gasteranthus quitensis | Oct | 22.2 | 12.1 | 12.9 | 12.3 | 12.8 | 12.0 | 12.5 | 12.1 |
| Gasteranthus quitensis | Sep | 20.0 | 12.2 | 11.9 | 12.2 | 11.5 | 11.6 | 11.9 | 13.1 |
| Glossoloma oblongicalyx | Apr | 0.0 | 18.7 | 17.3 | 17.0 | 17.5 | 17.1 | 17.5 | 16.8 |
| Glossoloma oblongicalyx | Aug | 35.3 | 18.6 | 21.9 | 22.6 | 21.4 | 22.7 | 21.2 | 23.5 |
| Glossoloma oblongicalyx | Dec | 0.0 | 18.0 | 18.2 | 18.0 | 17.8 | 18.3 | 18.9 | 18.4 |
| Glossoloma oblongicalyx | Feb | 0.0 | 18.1 | 18.6 | 20.2 | 19.5 | 19.4 | 18.9 | 19.4 |
| Glossoloma oblongicalyx | Jan | 0.0 | 17.4 | 14.4 | 13.3 | 14.6 | 14.0 | 14.4 | 13.4 |
| Glossoloma oblongicalyx | Jul | 38.9 | 18.6 | 18.9 | 21.0 | 19.6 | 19.7 | 20.1 | 20.7 |
| Glossoloma oblongicalyx | Jun | 10.0 | 18.4 | 19.8 | 20.2 | 20.3 | 20.2 | 19.6 | 20.1 |
| Glossoloma oblongicalyx | Mar | 0.0 | 19.0 | 14.6 | 14.0 | 13.7 | 14.7 | 14.7 | 13.6 |
| Glossoloma oblongicalyx | May | 11.8 | 18.4 | 15.4 | 14.1 | 15.7 | 14.6 | 16.5 | 14.8 |
| Glossoloma oblongicalyx | Nov | 22.2 | 18.4 | 20.5 | 21.0 | 20.2 | 22.9 | 19.7 | 23.2 |
| Glossoloma oblongicalyx | Oct | 38.9 | 18.5 | 21.1 | 21.4 | 20.5 | 21.9 | 20.4 | 22.9 |
| Glossoloma oblongicalyx | Sep | 35.0 | 18.7 | 17.8 | 18.6 | 18.1 | 17.6 | 17.9 | 18.0 |
| Glossoloma purpureum | Apr | 6.7 | 13.4 | 14.8 | 15.7 | 15.4 | 14.4 | 15.5 | 15.4 |
| Glossoloma purpureum | Aug | 17.6 | 13.9 | 10.3 | 10.6 | 10.8 | 10.5 | 12.1 | 10.0 |
| Glossoloma purpureum | Dec | 12.5 | 13.1 | 13.2 | 13.4 | 12.6 | 13.9 | 13.4 | 12.5 |
| Glossoloma purpureum | Feb | 12.5 | 13.2 | 13.7 | 13.7 | 13.8 | 13.5 | 14.0 | 13.5 |
| Glossoloma purpureum | Jan | 16.7 | 14.3 | 14.1 | 14.1 | 13.2 | 15.2 | 13.9 | 13.9 |
| Glossoloma purpureum | Jul | 22.2 | 13.7 | 14.8 | 15.4 | 15.2 | 14.9 | 15.9 | 14.2 |
| Glossoloma purpureum | Jun | 10.0 | 12.9 | 13.8 | 14.0 | 14.0 | 13.8 | 14.3 | 14.0 |
| Glossoloma purpureum | Mar | 0.0 | 13.4 | 13.9 | 13.9 | 14.0 | 14.8 | 13.3 | 14.3 |
| Glossoloma purpureum | May | 5.9 | 12.7 | 16.3 | 18.7 | 16.8 | 17.8 | 16.9 | 17.9 |
| Glossoloma purpureum | Nov | 11.1 | 14.6 | 10.7 | 8.9 | 10.0 | 8.8 | 10.8 | 8.6 |
| Glossoloma purpureum | Oct | 11.1 | 13.2 | 9.9 | 8.5 | 9.3 | 8.7 | 10.4 | 8.3 |
| Glossoloma purpureum | Sep | 25.0 | 13.9 | 15.2 | 15.6 | 15.2 | 15.2 | 15.0 | 15.3 |
| Kohleria affinis | Apr | 46.7 | 15.5 | 15.9 | 15.1 | 15.1 | 16.4 | 15.6 | 15.8 |
| Kohleria affinis | Aug | 5.9 | 15.0 | 16.2 | 16.6 | 17.0 | 16.2 | 17.1 | 15.7 |
| Kohleria affinis | Dec | 12.5 | 15.8 | 17.6 | 17.3 | 17.6 | 16.4 | 17.8 | 16.9 |
| Kohleria affinis | Feb | 25.0 | 15.4 | 15.0 | 14.4 | 15.6 | 14.8 | 15.1 | 15.4 |
| Kohleria affinis | Jan | 8.3 | 14.9 | 16.6 | 16.7 | 15.2 | 15.7 | 16.6 | 16.1 |
| Kohleria affinis | Jul | 0.0 | 15.4 | 14.4 | 14.0 | 13.8 | 15.5 | 14.3 | 14.2 |
| Kohleria affinis | Jun | 15.0 | 15.3 | 14.7 | 14.6 | 15.8 | 15.8 | 15.2 | 15.6 |
| Kohleria affinis | Mar | 41.7 | 16.1 | 16.1 | 16.0 | 16.5 | 15.6 | 15.8 | 15.6 |
| Kohleria affinis | May | 29.4 | 15.7 | 15.5 | 16.3 | 17.2 | 16.8 | 16.5 | 15.4 |
| Kohleria affinis | Nov | 0.0 | 15.4 | 16.8 | 15.6 | 17.5 | 15.5 | 17.2 | 17.2 |
| Kohleria affinis | Oct | 0.0 | 14.9 | 17.1 | 16.6 | 17.3 | 16.2 | 16.8 | 16.4 |
| Kohleria affinis | Sep | 5.0 | 15.1 | 15.1 | 15.5 | 15.4 | 14.9 | 15.6 | 15.8 |
Grab an example Drymonia teuscheri v Drymonia collegarum
D["Drymonia teuscheri","Drymonia collegarum"]
## [1] 0.02974
Dtraits["Drymonia teuscheri","Drymonia collegarum"]
## [1] 0.04191869
Dint["Drymonia teuscheri","Drymonia collegarum"]
## [1] 0.4755805
Logit E
inv.logit E
## # A tibble: 6 x 4
## # Groups: Model [6]
## Model Var1 Var2 Correlation_E
## <chr> <fct> <fct> <dbl>
## 1 interaction_attraction Drymonia teuscheri Drymonia colle… -0.492
## 2 interaction_repulsion Drymonia teuscheri Drymonia colle… 0.0870
## 3 phylogenetic_attraction Drymonia teuscheri Drymonia colle… 0.784
## 4 phylogenetic_repulsion Drymonia teuscheri Drymonia colle… -0.864
## 5 trait_attraction Drymonia teuscheri Drymonia colle… 0.698
## 6 trait_repulsion Drymonia teuscheri Drymonia colle… -0.867
## # A tibble: 7 x 4
## # Groups: Model [7]
## Model Var1 Var2 Correlation_P
## <chr> <fct> <fct> <dbl>
## 1 baseline Drymonia teuscheri Drymonia colle… NA
## 2 interaction_attraction Drymonia teuscheri Drymonia colle… -0.394
## 3 interaction_repulsion Drymonia teuscheri Drymonia colle… -0.160
## 4 phylogenetic_attraction Drymonia teuscheri Drymonia colle… 0.495
## 5 phylogenetic_repulsion Drymonia teuscheri Drymonia colle… -0.469
## 6 trait_attraction Drymonia teuscheri Drymonia colle… 0.355
## 7 trait_repulsion Drymonia teuscheri Drymonia colle… -0.471
## # A tibble: 7 x 4
## # Groups: Model [7]
## Model Var1 Var2 Correlation_Ynew
## <chr> <fct> <fct> <dbl>
## 1 baseline Drymonia teuscheri Drymonia co… 0.0889
## 2 interaction_attraction Drymonia teuscheri Drymonia co… -0.394
## 3 interaction_repulsion Drymonia teuscheri Drymonia co… -0.162
## 4 phylogenetic_attraction Drymonia teuscheri Drymonia co… 0.490
## 5 phylogenetic_repulsion Drymonia teuscheri Drymonia co… -0.468
## 6 trait_attraction Drymonia teuscheri Drymonia co… 0.355
## 7 trait_repulsion Drymonia teuscheri Drymonia co… -0.471